GPU dynamic memory allocation - #18

#!/bin/bash
set -eu
printf '%s\n' '--- candidate files ---'
git ls-files | grep -E '(^|/)(DeviceAllocator\.h|sycl\.h)$' || true
printf '%s\n' '--- relevant symbols and call sites ---'
rg -n -C 6 'class Stack|struct Stack|void init|parallel_for|syclItem|get_nd_item' include/DeviceAllocator.h include/sycl.h 2>/dev/null || true
printf '%s\n' '--- file sizes ---'
wc -l include/DeviceAllocator.h include/sycl.h 2>/dev/null || true

#!/bin/bash
set -eu
printf '%s\n' '--- DeviceAllocator.h surrounding implementation ---'
cat -n include/DeviceAllocator.h | sed -n '1,150p'
printf '%s\n' '--- init call sites ---'
rg -n -C 5 '\.init\(\)|init\(\)' --glob '!include/DeviceAllocator.h' --glob '!include/sycl.h' .
printf '%s\n' '--- parallel_for signatures in repository ---'
rg -n -C 3 'parallel_for\s*\(' --glob '*.h' --glob '*.hpp' --glob '*.cpp' --glob '*.cc' --glob '*.cxx' .

#!/bin/bash
# Find device-side default construction of arrays that use GlobalDeviceAllocator.
rg -nP -C4 '\barray\s*<[^>]*(GlobalDeviceAllocator|Allocator<)[^>]*>\s+\w+\s*(;|\{\s*\})' --type=cpp --type=cc -g '!**/build/**'
# Inspect the terminating allocator and allocate() paths for null handling.
fd -t f 'DeviceAllocator.h' | xargs -r sed -n '80,200p'

#!/bin/bash
set -eu

echo '--- GlobalDeviceAllocator and array references ---'
rg -n -C3 'GlobalDeviceAllocator|m_allocator|array\s*\(' include src test examples --glob '!**/build/**' 2>/dev/null || true

echo '--- arrays.h structure and constructors ---'
fd -t f 'arrays.h' . | while read -r f; do
  echo "FILE: $f"
  ast-grep outline "$f" | head -120
  rg -n -C8 'm_allocator|array\s*\(|allocate\s*\(' "$f"
done

echo '--- array declarations and allocator aliases ---'
rg -n -C3 'using .*Allocator|typedef .*Allocator|array\s*<' . \
  --glob '*.{h,hpp,c,cpp,cc,cxx}' \
  --glob '!**/build/**' \
  --glob '!**/dist/**' 2>/dev/null || true

#!/bin/bash
set -eu

echo '--- All direct GlobalDeviceAllocator uses ---'
rg -n -C5 'GlobalDeviceAllocator' . --glob '!**/build/**' --glob '!**/dist/**'

echo '--- SYCL allocator aliases and device array declarations ---'
fd -t f -e h -e hpp -e cpp -e cc -e cxx . include src test models | while read -r f; do
  if rg -q 'DeviceAllocator|GlobalDeviceAllocator|__SYCL_DEVICE_ONLY__|array_ns::array|array<' "$f"; then
    echo "FILE: $f"
    rg -n -C4 'DeviceAllocator|GlobalDeviceAllocator|__SYCL_DEVICE_ONLY__|array_ns::array|array<' "$f" | head -160
  fi
done

echo '--- sycl.h relevant definitions ---'
fd -t f 'sycl.h' . | while read -r f; do
  echo "FILE: $f"
  sed -n '1,180p' "$f"
done

echo '--- arrays.h include and allocator-selection context ---'
sed -n '1,80p' include/arrays.h
sed -n '1240,1300p' include/arrays.h

#!/bin/bash
set -eu

echo '--- ModelData and EcolabPoint definitions ---'
sed -n '1,125p' models/ecolab_model.h

echo '--- DeviceType definition and uses ---'
rg -n -C8 'template.*DeviceType|struct DeviceType|class DeviceType|DeviceType<' include models src \
  --glob '*.{h,hpp,c,cpp,cc,cxx}' | head -240

echo '--- Device-side default construction of GlobalDeviceAllocator-backed objects ---'
rg -n -C8 'UnsignedArray|Allocator<|cell_ids|density|array<[^;]*(GlobalDeviceAllocator|ModelData::Allocator)' \
  models/ecolab_model.h models/ecolab_model.cc include \
  --glob '*.{h,hpp,c,cpp,cc,cxx}' | head -320

echo '--- Read-only semantic probe of the relevant preprocessor branches ---'
python3 - <<'PY'
from pathlib import Path

allocator = Path("include/DeviceAllocator.h").read_text()
arrays = Path("include/arrays.h").read_text()
model = Path("models/ecolab_model.h").read_text()

checks = {
    "device constructor initializes allocator to nullptr":
        "GlobalDeviceAllocator(): allocator(nullptr) {}" in allocator,
    "device allocate dereferences allocator without a guard":
        "{return reinterpret_cast<T*>(allocator->allocate(n*sizeof(T)));}" in allocator,
    "array default constructor is defaulted":
        "array()=default;" in arrays,
    "array stores allocator as a value member":
        "A m_allocator;" in arrays,
    "EcolabPoint uses GlobalDeviceAllocator under SYCL":
        "template <class T> using Allocator=GlobalDeviceAllocator<T>;" in model,
    "EcolabPoint has a GlobalDeviceAllocator-backed array member":
        "array<int,Allocator<int>> density;" in model,
}
for name, result in checks.items():
    print(f"{'PASS' if result else 'FAIL'}: {name}")
PY

#!/bin/bash
set -euo pipefail

# Map nsp declarations and definitions before inspecting call sites.
ast-grep outline models --items all --type function --match nsp

# Inspect C++ call sites and nearby const-qualified model references.
rg -n --type=cpp -C4 '(?:\.|->)nsp\s*\(|const\s+SpatialModel' .

#!/bin/bash
set -euo pipefail

printf '%s\n' '--- declarations and definitions ---'
rg -n -C3 '\bnsp\s*\(' models --glob '*.{h,hpp,hh,cc,cpp,cxx}' --glob '!*.md'

printf '%s\n' '--- all nsp call sites ---'
rg -n -C3 '\bnsp\s*\(' . --glob '*.{h,hpp,hh,cc,cpp,cxx,c}'

printf '%s\n' '--- SpatialModel type uses ---'
rg -n -C3 '\b(const\s+)?SpatialModel\b' . --glob '*.{h,hpp,hh,cc,cpp,cxx}'